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I. INTRODUCTION 


In recent years, several new advancements have been 
made in the area o£ navigation to enable the pilot to 
navigate more accurately and with greater efficiency. 

These advancements have arisen because of the require- 
ment for improved enroute navigation and landing systems 
to provide increased safety and greater capability for 
instrument flight rule (IFR) flights in adverse weather 
conditions. The desire for flight capability during IFR 
conditions arises from increased demand for air travel 
under all weather conditions. Furthermore, increased 
air traffic densities place greater demands on safety 
which has a potential impact on improved navigation 
accuracy requirements. 

The equipment modifications to provide the navig- 
ation improvements range from airframe avionics changes 
for implementation of new enroute navigation techniques, 
to modified landing systems which provide greater pre- 
cision in final approach and landing. Two such systems 
that this report is concerned with are the area navig- 
ation (RNAV) concept and the Microwave Landing System 
(MLS) . Both of these concepts are elements of the 
National Airspace Systems (NAS) Upgraded Third Generation 
Plan. The area navigation system concept allows for 
flight between predetermined navigation waypoints not 
necessarily coincident with VOR/DME ground stations util- 
ized if, in fact, the RNAV support is supplied by the 
VOR/DME navigation system. RNAV can also be supported 
by other navigation systems such as LORAN-C, Omega, and 
GPS which are more amenable to the point-to-point navig- 
tion concept associated with RNAV. The Microwave Landing 
System is an improved landing aid which can provide 4-D 
guidance with accuracies sufficient to provide Category 
III-C landings. With suitable avionics the MLS mode of 
navigation can be utilized as a highly accurate RNAV 
system in the terminal area. 

To date, many studies have been performed that con- 
centrate on either the area navigation concept or the 
Microwave Landing System. However, very little has been 
done in the way of examining the aircraft transition 
from the RNAV region to the MLS region. Although both 
systems are capable of providing three-dimensional position 
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information, the accuracies of the two systems are quite 
different. These differences may produce considerable 
transients, or the RNAV delivery accuracies at the MLS 
coverage boundary may be unacceptably large, under cer- 
tain conditions. These potential problems provide the 
motivation for the study documented in this report. 

Another motivating factor is to provide analytical 
support to the on-going terminal configured vehicle (TCV) 
program being conducted at the NASA Langley Research Center. 
This study includes evaluation of the RNAV/MLS transition 
problems. The TCV aircraft is a Boeing 737 equipped with 
an inertial RNAV system updated with VOR/DME or DME/DME. 
Sensors and software are also available onboard to pro- 
vide MLS guidance. 

These motivating factors, form the basis in establishing 
the objectives of this avionics sensitivity study. The 
primary objective is to generate statistical data re- 
garding aircraft position errors in the RNAV/MLS transition 
region. This data is to be generated for the TCV aircraft 
with the navigation components aboard this particular air- 
craft. The statistical data is to be generated such that 
the sensitivity to navigation errors and flight geometries 
can be determined. 

The primary purpose for generating the data is to 
provide insight into the behavior of the aircraft in the 
RNAV/MLS transition region. Knowledge of aircraft behavior 
in this region due to navigation errors is desirable for 
subsequent guidance law and airspace design. The avionics 
sensitivity data will provide the information necessary to 
establish requirements for additional control law design 
during transition and to establish airspace requirements 
for maneuvering to null out any residual RNAV errors upon 
MLS transition. Availability of the statistical data 
will also provide insight into the evaluation of the feas- 
ibility of RNAV/MLS transition for various geometries, 
flight profiles and error parameter magnitudes. Feasibility 
of the various profiles and geometries aid in the establish- 
ment of guidelines for route design in the RNAV/MLS trans- 
ition region. Finally, the data base is beneficial for 
providing insight into the development of subsequent flight 
test plans. 

To meet the stated objectives, a computer simulation 
developed by Systems Control, Inc. (Vt) [SCI (Vt) ] was 
appropriately modified. The digital simulation provides 
data regarding the expected aircraft position errors along 
any route or route position for the TCV aircraft in both the 
RNAV and MLS regions. It is possible, with the digital 
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simulation, to evaluate navigation system performance for 
any operational scenario and range of error parameters. In 
this manner, the feasibility of proposed navigation and 
facility locations and approach path/runway configurations 
can be evaluated in a cost-effective manner. 

The flight profiles of the current evaluation of navi- 
gation system performance are restricted to flight paths 
typical of conventional take-off and landing (CTOL) air- 
craft. However, navigation facility locations are selected 
to give worst-case error possibilities. The error sources 
have been selected to be compatible with the recommended 
error budget values as discussed in AC 90-45A [1]. The 
relative amount of random bias and noise contained within 
the error budget values, however, is varied so that it is 
possible to evaluate the effect of each. 

Through judicious quantification of the variables 
in the study, a data base is generated for subsequent 
applications and analysis. The data is basically pre- 
sented in the form of plots and graphs. The organization 
of the report is as follows. The second section presents 
the motivation for implementing an ensemble or covariance 
type digital simulation as opposed to a Monte Carlo type. 
This discussion includes a presentation of the advantages 
and disadvantages of each of the simulation techniques. 

The discussion leads into the presentation of the RNAV and 
MLS simulations. Supporting analysis is presented in 
Appendix A. Chapter II ends with a discussion of the 
simulation capabilities. 

The analytical approach is presented in Chapter III. 

The section begins with a discussion of the application of 
the ensemble simulation to meet the study objectives as 
stated in this chapter. To meet the study objectives it 
is necessary to carefully design a suitable experimentation 
test plan. The test plan for the study documented in this 
report is presented in the remainder of Chapter III with 
supporting details presented in Appendix B. 

The analytical results of the sensitivity study are 
presented in Chapter IV. The primary quantities of interest 
are the cross-track and vertical deviation time histories 
through the RNAV/MLS transition as a function of the various 
study parameters. These time histories are presented in 
the form of plots with corresponding narrative discussing 
the contents of each. 

The study conclusions are presented in Chapter V with 
corresponding recommendations presented in Chapter VI . 
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Generally, two techniques are associated with simul- 
ation experimentation. These are Monte Carlo analysis 
and ensemble analysis (covariance propagation). The former 
represents the occurrence of a single random event where- 
as the latter represents an ensemble of events occurring 
simultaneously. For example, in a Monte Carlo simulation, 
a single aircraft is executing a simulated flight along 
a prescribed flight path driven by errors representative 
of a single source. In the ensemble simulation, an en- 
semble of aircraft are traversing the prescribed flight 
path driven by errors representing the ensemble average 
of a particular source. The primary difference between 
these two techniques is that the analytical models 
of the Monte Carlo simulation can be either linear or 
nonlinear, whereas the analytical models of the ensemble 
simulation must be linearized. Hence, the Monte Carlo 
simulation can be utilized to evaluate all effects whether 
linear or nonlinear. Judicious modelling must be employed 
for the linearized ensemble simulation to retain realism. 
The primary advantage of the ensemble simulation is that 
statistical data is achievable through a single simulation 
execution, whereas, similar statistical data is achievable 
from the Monte Carlo simulation only through a large number 
of computer runs. The level of confidence in the Monte 
Carlo generated statistics is directly dependent on the 
sample size. Hence, this approach is potentially more 
costly than using covariance propagation techniques, es- 
pecially for the generation of preliminary statistical 
data. 


For this study the linearized ensemble approach was 
utilized to generate the desired statistical data. For 
the system under investigation the linearization assump- 
tion is not overly restrictive. This is due to the fact 
that suitable linear navigation error models exist and 
the navigation filters and the controls already exist 
in linear form. Furthermore, since this is essentially 
a navigation system evaluation, the aircraft dynamics 
and autopilots can be represented by simplified linear 
models. The primary nonlinearities omitted in the analysis 
are the control bounds. However, these nonlinearities 
will not appreciably affect the ensemble results. 

The sensitivity analysis of this study is concerned 
with aircraft performance during transition from an 
area navigation (RNAV) environment to the Microwave Landing 
System (MLS) coverage region. The ensemble simulation 
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development therefore involves a separate analysis for 
the RNAV and MLS avionics. Transition is achieved by 
retaining those elements of the covariance matrix common 
to both systems (aircraft dynamics, autopilot, and so on) 
at MLS lock-on and substituting the RNAV-related elements 
with MLS-related elements. In this manner the covariance 
propagation across the MLS coverage boundary will represent 
the actual situation. This approach precludes any potential 
blending of RNAV and MLS signals. Therefore the discussion 
regarding simulation development is presented in two 
parts - RNAV and MLS simulation development. 

Generally, the analytical models were derived from 
existing literature on the TCV-737 aircraft system [1-4]. 
Where insufficient information was available to derive 
adequate analytical models other suitable sources were 
consulted [5-7] (INS, VOR, DME, MLS). As indicated prev- 
iously, almost all of the navigation filter and the control 
equations (except for control bounds) appear in linear 
form, hence, simplifying the linearization process. The 
aircraft dynamics are assumed to represent simple point 
mass motion and are represented by second order equations 
of motion in both the vertical and lateral directions. 

The section organization is as follows. Section 2.2 
presents a detailed discussion of the simulation model 
development. This includes the models of the RNAV and 
MLS portions of the simulation. Block diagrams are pre- 
sented to indicate information flow in these two systems. 

The assumptions and derivation of the covariance propa- 
gation equations are then presented. Section 2.3 des- 
cribes in detail the features of the simulation design. 

These features include run time waypoint and station speci- 
fication, interactive control logic, free format input, 
printer and CALCOMP plotting capabilities and others. 
Finally, Section 2.4 presents a summary of the overall 
simulation capabilities. 

2.2 SIMULATION MODELS 

2.2.1 RNAV Simulation 


As indicated in Section 2.1, the basic TCV-737 linear 
ensemble simulation consists of two parts. The first is 
the RNAV simulation and the second is the MLS simulation. 
The purpose of this subsection is to present, with the aid 
of a block diagram, the details of the RNAV simulation. 

The details with regard to mathematical expressions for 
each component of the simulation are presented in Appendix 
A. 
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Figure 2.1 is the block diagram describing the RNAV 
simulation. The components o£ this particular simulation 
include an Inertial Navigation System (INS) , second order 
lateral and vertical complementary filters, lateral and 
vertical autopilots, barometric altimeter and aircraft 
dynamics. The lateral inertial velocity, and 

6Vgi, and vertical inertial acceleration, , serve 

as inputs which complement the second order lateral and 
vertical filters respectively. The filter estimated 
north and east position errors, and , are 

transformed to station-referenced range, Ap^ , and 

bearing, A9^ , to yield the computed measurement errors. 

The geometry depicting the range and bearing measurements 
is shown in Figure 2.2. The computed measurement errors 
are subtracted from the actual measurement errors and this 
differenced quantity is filtered after a coordinate 
transformation into the second order lateral complementary 
filter. The filter estimated lateral position and velocity 
states are u§ed to compute the indjcated^cross- track 
deviation CTD and heading error, \Ij . CTD and $ are 
used to compute the roll command ^ which is input to 
the aircraft dynamics. 

In the vertical channel the input to the complementary 
filter is the measured altitude and altitude rate as 
output from the barometric altimeter model. This model 
includes bias and noise input to a first order lag to 
represent the time lag characteristic of altimeters. The 
output is a commanded altitude rate which, together with 
the measured altitude form the basic input into the vertical 
autopilot. For 3D RNAV systems the along- track error prop- 
agates into the vertical channel as shown in Figure 2.3. 

The prim=xry impact of the along-;, rack error is to shift the 
vertical profile. This error source can be significant 
for 3D RNAV systems, especially when executing steep 
descent profiles (for example, a 6® STOL approach). This 
is due to the fact that the vertical error is proportional 
to the glideslope according to 

VD = ATD tan y C^.l) 

Hence, for an along-track error of 2780 meters (which is 
reasonable for RNAV systems) and a glideslope of 3°, the 
vertical deviation due to along-track alone is 146 meters. 
The output of the vertical autopilot is commanded vertical 
acceleration (pitch rate) which is fed into the vertical 
dynamic equations. 
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Figure 2.1 VOR/DME RNAV Simulation Block Diagram 
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Figure 2.3 Along-Track Deviation Contribution to 
Total Vertical Deviation 


For the current simulation the lateral autopilot and 
dynamics and the vertical autopilot and dynamics are each 
represented by second-order equations. The aircraft 
dynamic output is the actual cross -track and vertical devia- 
tion which is fed back to evaluate the measurements. 

The error sources for the RNAV simulation consist 
of range and bearing noise and bias, altimeter noise and 
bias and INS error (platform tilt errors, accelerometer 
biases and drifts) . For this effort the INS and altimeter 
errors will be fixed. 

2.2.2 MLS Simulation 


The flow chart depicting the MLS simulation is shown 
in Figure 2.4. The lateral and vertical autopilot and 
the aircraft dynamics are essentially the same as those 
of the RNAV simulation. The primary difference in the 
two simulations is the manner in which the guidance comput- 
ations are computed from which appropriate steering commands 
can be generated. 

The guidance computations are based on MLS-derived 
position measurements. The MLS measurements differ from 
the RNAV measurements in that they consist of range, azi- 
muth angle and elevation angle measured with respect to 
a runway fixed coordinate system. Figure 2.5 schematically 
depicts the MLS coverage region and Figure 2.6 the corres- 
ponding measurements. Hence, several coordinate transfor- 
mations are required to properly perform the guidance 
computations . 


11 



^S> 



Figure 2.4 MLS Simulation Block Diagram 
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MLS Coverage Region 






In the MLS mode, the second order complementary filter 
has as input the three acceleration components. The 
accelerations are derived from body fixed accelerometers. 
Runway- referenced accelerations are derived using aircraft 
attitude information derived from the attitude/heading 
reference unit to transform the accelerometer acceleration 
components. The filtered estimated position errors, 

Xr, yj^, and (computed measurements) are subtracted 

from the MLS- derived measurements to form the second input 
to the second-order filter. 

The filter estimated velocities are used to compute 
heading errors. Similarly, the filtered estimated posi- 
tions are used to compute course deviations. The heading 
error and course deviations are then used to generate the 
appropriate steering commands. 

The error sources for the MLS mode of operation con- 
sist of the measurement errors and errors associated with 
the acceleration measurements. The MLS measurement errors 
consist of range, azimuth and elevation bias and noise. 
Similarly, the alignment and accelerometer errors are repre- 
sented by bias and random noise. 

Currently, blending between the RNAV and MLS measure- 
ments is not envisioned. Hence, at the MLS boundary 
the navigation mode is assumed to transfer discretely 
from RNAV to MLS. The reason for this is that not much 
is known about blending techniques and it is uncertain 
what effects are to be expected at MLS transition from 
RNAV. The purpose of this study is to provide more insight 
into this problem area so that subsequent analysis can 
be performed with regards to the development of approp- 
riate blending techniques. 

2.2.3 Covariance Propagation 

The previous two subsections have discussed the two 
basic simulations that make up the overall ensemble. - 
simulation for this study. Appendix A has presented the 
equations describing the various components of the RNAV 
and MLS simulations. These equations are presented in 
linearized form suitable for ensemble simulation implement- 
ation. The linearity assumption is the most important 
one regarding covariance propagation yet it is this very 
assumption that limits the applicability of ensemble analy- 
sis. This assumption eliminates the evaluation of system 
performance due to inherent nonlinearities unless these 
nonlinearities can be linearized without destroying real- 
ism or if they can be represented by such techniques as 


15 



describing functions. However, ensemble techniques do 
provide the capability of generating desirable statistical 
data in a cost effective manner. It is this cost effective- 
ness that provides the motivation for utilizing this 
approach. 

For a dead-reckoning navig'ation system, such as 
an inertial or air data system, updated with external 
measurements, VOR/DME or DME/DME, the covariance propa- 
gation can be performed in one of two ways. First of all 
the covariance matrix can be propagated between ground- 
based measurement updates using the dynamics associated 
with the dead- reckoning system, autopilot and aircraft. 

The covariance would then be updated at the measurement 
interval using the appropriate measurement updating tech- 
niques. The covariance propagation for this approach 
is represented by (in discrete form) [8] 


State: x^ = (|)j^ 

-1 ^k-i 

+ W^_l, N(0,Q^) 

(2.2) 

Measurement: z 

k ' »k 

Xk + Vj^'v. N(0,Rj^) 

(2.3) 

State Estimate 




Extrapolation : 

X^(-) 

" ‘^k-l ^k-1^'*'^ 

(2.4) 

Covariance „ 

Extrapolation: 

Pkt-’ 

“ ‘•’k-1 ^k-l^"^^ k-1 ^k- 

l(2.5) 

State Estimate 

Update : 



(2.6) 

Covariance „ 

Update : 

PkCO 

= [i-WPkt-ni-W‘ 

* wl 

+ 

(2.7) 


where : 

X - actual state 

(|) - state transition matrix 

W - process noise 
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The states, x, state transition matrix, (|>, and 
process noise matrix, Q, are presented explicitly in 
Appendix A. 

2.3 BASIC SIMULATION DESIGN 

The simulation developed for this study was intention- 
ally made comprehensive and flexible. It is comprehensive 
in that complete guidance, navigation, and control systems 
are included, along with models of the full complement of 
radio and inertial navigation sensors. Flexibility is 
achieved through modularity of design in conjunction with 
a sophisticated executive structure and user interface, 
which had previously been developed and verified. The 
executive structure provides for employing multiple iter- 
ation rates of the simulator functions Cguidance, navigation, 
control, covariance propagation, output recording, etc.) 
and for accommodating large simulator core requirements in 
an efficient manner. The user interface permits free- 
format input/output. The simulation also has an interactive 
feature whereby the engineer may monitor, stop, correct, 
and resume simulation computations from a remote terminal; 
this feature was used extensively in the simulation develop- 
ment process as an efficient debugging aid. 

The simulation has been designed to operate in two 
modes. The first is a nonlinear Monte Carlo mode which 
provides the capability for evaluation of nonlinear effects 
on a case-by-case basis. The second is the covariance prop- 
agation (ensemble) mode which provides the capability of 
evaluating the ensemble statistics on a single computer run 
basis. Execution of a large number of Monte Carlo runs 
will approach the statistics generated by the ensemble simul- 
ation for the same error source inputs. 

2.4 SIMULATION CAPABILITIES 

The purpose of this subsection is to summarize the 
ensemble simulation capabilities. As indicated in the 
introduction the primary objective for this study is to 
propagate aircraft position errors in the RNAV/MLS trans- 
ition region. This objective is the primary driving 
force in the determination of the specific simulation 
capabilities listed here. 

Although the primary application for this study is 
the ensemble mode it is necessary to generate the nominal 
trajectory. For this study the nominal is generated as 
an error free Monte Carlo simulation. Hence, by adding 
errors it is possible to determine the mean three-dimensional 
position errors from constant or repetitive error sources. 
This part of the simulation includes such nonlinearities as 
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z - measurements 
H - measurement matrix 
V - measurement noise 
X - state estimate 
P - error covariance matrix 

Q - diagonal matrix representing process noise 
K - weighting matrix 
I - identity matrix 

R - diagonal matrix representing measurement noise 

For this particular study the errors associated with 
the inertial navigation system and altimeter would appear 
as process noise. The measurements would be VOR/DME, DME/ 
DME or MLS. The errors associated with each o£ the measure- 
ments would be represented as measurement noise, for the 
random component, and as additional states for the bias 
or correlated noise component. 

In the TCV-737 navigation system the inertial navi- 
gation system is updated every second with external pos- 
ition fixes. In the interest of conserving computation 
time it is desirable to integrate with a step size of 
at least one second or, if the system time constants per- 
mit, at a larger step size. This can be accomplished by 
including the measurements as part of the state and their 
associated noise components as process noise. Hence, the 
covariance extrapolation is reduced to 

‘^k-1 ^k-1 '*‘Ll ^ ^k-1 ^ 

where includes inertial navigation system, alti- 

meter, VOR, DME, or MLS noise components as applicable. 

Note that the block diagrams for the RNAV and MLS simul- 
ations of the previous sections have been configured in 
a form suitable for application of the above covariance 
extrapolation. Utilizing this particular formulation 
negates the capability of examining the effect of update 
rate. However, for the current study this disadvantage 
is insignificant compared to the advantage gained by 
reduced computation time. 
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control bounds. The ensemble portion of the simulation 
provides the capability of generating the standard devia- 
tion of the position errors due to random errors about the 
mean. The ensemble simulation, due to the inherent linear- 
ity assumptions in the development of this type of simula- 
tion, provides the capability of evaluating the statistical 
behavior for a larger range of parameters through application 
of the root-sum-squaring (RSS) technique to the simulation- 
generated results. 

The error propagation capabilities provide the 
analytical tool necessary to determine the sensitivity 
of position errors to variations in several key parameters. 
The first parameter of interest is the error source type 
and magnitude. The type can be, for example, bias, cor- 
related or random noise and the magnitude can represent 
current navigation accuracy requirements or proposed re- 
quirements . In this manner, hypothetical navigation 
system performance can be evaluated by appropriately 
varying the simulation inputs. 

The next parameter of interest, especially for angle 
measurement or range -dependent navigation systems, is the 
flight configuration (speed and geometry) . The simulation 
has the capability to have as input and prespecified three- 
dimensional route with arbitrary ground facility locations. 
This provides the ability to evaluate navigation system 
performance for any operational scenario. This flexibility 
is important for the avionics sensitivity study since the 
design of RNAV/MLS approach paths and the allocation of 
airspace depend to a great extent on the potential RNAV 
delivery errors at the MLS boundary and the capability to 
deal with them upon MLS transition. These errors are 
geometry-dependent and therefore it becomes necessary to 
evaluate the impact of geometry on the achievable position 
accuracies . 

The capabilities of the navigation simulation can be 
summarized as; 

(1) determination of the mean three-dimensional posi- 
tion error from constant or repetitive error 
sources ; 

(2) determination of the standard deviation of the 
position errors due to random error sources about 
the mean; 

(3) combination of (1) and (2) to obtain the "mean 
square error"; 
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(4) error distribution using RSS technique; 

(5) sensitivity o£ position errors to variations in 
error source types and magnitudes; 

(6) sensitivity of position errors to variations in 
flight configuration (speed and geometry); 

(7) response of the system to nonlinear error sources; 

(8) comparison of filters for error attenuation; 

(9) comparison of RNAV and MLS system mechanizations; 

(10) performance of the system using real navigation 
aid data; 

(11) determination of RNAV and MLS system performance 
along any prespecified route; and 

(12) propagation of non-Gaussian error probability 
distributions through typical RNAV and MLS systems. 



III. ANALYTICAL APPROACH 


3.1 INTRODUCTION 

The basis of an avionics sensitivity study is the 
generation of a suitable statistical data base. The 
primary objective is to evaluate the propagation of 
navigation error statistics through the system and 
assess the impact on aircraft position performance. 

Several analytical approaches can be taken to generate 
the desired data base. The first approach is using 
Monte Carlo type simulations and the second is the use of 
covariance propagation (ensemble simulation) . The dif- 
ference between these two simulation techniques has been 
discussed in Section 2.1. The purpose of the discussion 
in this section is to indicate the manner in which the 
ensemble simulation is applied to the current avionics 
sensitivity study. 

The ensemble simulation (covariance propagation) 
essentially generates the time history of the variances 
and covariances of the states along a prespecified nominal 
trajectory. The covariance time history generated in 
this manner eliminates the necessity to generate a large 
number of Monte Carlo runs. The limitation, as discussed 
previously, is the requirement to linearize the approp- 
riate system equations. For this study the primary non- 
linearities exist in the equations of motion for the air- 
craft and in the control bounds. However, the focus here 
is on avionics sensitivities to navigation errors, hence, 
the aircraft dynamics play a minor role and can be simpli- 
fied to linear form. 

The primary objective for implementation of the en- 
semble simulation is to generate as broad a statistical 
data base as feasible in the most cost effective manner. 
Through judicious selection of states and the equations 
governing their motion this objective can be achieved. 
Where the equations are not in linear form, a lineariz- 
ation process is executed such that the appropriate co- 
variance propagation form is achieved while retaining 
adequate physical representation. 

Following the development of the appropriate covari- 
ance propagation equations, it is necessary to generate a 
simulation test plan designed to include a sufficiently 
broad scope of parameters and broad range of variation of 
these parameters. The purpose of the data generated is to 
provide increased insight into the transient behavior of 
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navigation- induced errors in the RNAV/MLS transition 
region. Therefore, the test plan should include at least 
potentially worst case geometries and error parameter mag- 
nitudes and other geometries providing data of general 
interest. The test plan cannot be all encompassing, hence, 
the parameters and geometries for the sensitivity analysis 
must be judiciously selected. 

Application of the ensemble simulation developed 
for the avionics sensitivity analysis will yield expected 
aircraft position error time histories in the terminal 
area where navigation support is provided by RNAV and MLS. 

Of special interest is the expected transient at the MLS 
boundary due to the known error differences between the 
two navigation systems. The error differences are observed 
by comparing Tables 3.1 and 3.2. Table 3.1 indicates the 
MLS error specifications and Table 3.2 shows the VOR/DME 
errors. The VOR/DME system is currently the basic means 
of navigation support for RNAV. Other systems, such as 
LORAN-C, Omega and GPS, for example, are also typical RNAV 
systems. However, evaluation of aircraft position error 
due to these systems is beyond the scope of the current 
study . 

The primary purpose for the data base generated under 
this study is to provide appropriate information for guid- 
ance law development and airspace design. More specifically 
the data will, first of all, aid in establishing the feas- 
ibility of RNAV/MLS transition for various geometries, flight 
profiles and error parameter magnitudes. Having estab- 
lished the feasibility of RNAV/MLS transition the data 
provides information to establish requirements for addit- 
ional control laws during transition and to establish 
airspace requirements for maneuvering to null out RNAV 
errors. The airspace requirements would also include 
guidelines for route design in the RNAV/MLS transition. 
Finally, the data base can be used to identify potential 
areas or geometries for additional flight experimentation. 

Generation of a suitable data base for the avionics 
sensitivity analysis requires at the onset the development 
of a suitable simulation test plan. Because of the in- 
herent flexibility available for simulation experimentation, 
a broader range of parameters can be varied. For this 
study, these parameters are: flight profile, ground navi- 

gation aid location, runway (MLS) orientation, and error 
parameter content and magnitudes. These parameters and the 
range of values associated with each are selected to give 
as broad a data base as is feasible. They have been 
selected to yield results providing the greatest amount of 
information for subsequent analysis. 
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Table 3.1 


RTCA [9] MLS Specifications (le) ) 


Configuration 
Operational Use 

D 

Cat. I 

F 

Cat. II 

K 

Cat. Ill 

DME 

Bias 

91.4 m 

30.5 m 

6.1 m 


(300 ft.) 

(100 ft.) 

(20 ft.) 

Random 

* 

'k 


Total 

91.4 m 

30. 5 m 

6.1 m 

KL 




Bias 

.125 degrees 

.090 degrees 

.036 degrees 

Random 

.065 degrees 

.033 degrees 

.024 degrees 

Total 

.141 degrees 

.096 degrees 

.042 degrees 

EL 

Bias ' 

.050 degrees 

.050 degrees 

.050 degrees 

Random 

.058 degrees 

.035 degrees 

.035 degrees 

Total 

.077 degrees 

.061 degrees 

.061 degrees 


* Random error negligible compared to bias 


Table 3.2 

AC 90-4SA [10] VOR/DME Specifications (2o) 


SOURCE 

ERROR MAGNITUDE 

VOR; 


Ground 

1.9° 

Airborne 

O 

O 

CO 

DME: 


Ground 

185m (0.1 nmi ) 

Airborne 

926m (0.5 nmi) 


n 

















3.2 SIMULATION IMPLEMENTATION 


The purpose of this subsection is to explain in 
greater detail the implementation of the simulation and 
the issues alluded to in the introduction. One of the 
advantages of computer simulation is the capability to 
begin and end a particular simulated flight at arbitrary 
points. Since this study is primarily concerned with RNAV/ 
MLS transition and the evaluation of the expected trans- 
ients during this transition, it is desirable to initiate 
the simulation some time prior to the MLS boundary and 
terminate the simulation after the desired transients 
have been observed. The simulation must be initiated, 
however, with sufficient time prior to the MLS boundary 
to eliminate initial condition transients. Also, the 
termination of the simulation, for this study, is con- 
sidered to occur at the initiation of flare in order to 
observe the behavior of MLS-induced position errors during 
the entire flight in this navigation mode. 

Specific flight profiles and ground navigation facility 
locations are selected to provide sensitivity data valu- 
able for analysis and inputs for subsequent flight test 
design. The computer simulation flexibility is again util- 
ized for these parameter variations. Since the flight pro- 
file and navigation aid locations are varied by input 
into the simulation they need not be fixed as they are 
for actual flight tests. In this manner it is possible 
to design the flight profile, locate the navigation facil- 
ities and orient the runway and MLS antennas in such a 
manner to produce the most useful results. 

An additional degree of flexibility provided by the 
simulation approach is the capability to vary the error 
source content and magnitude. Variation of the error para- 
meters yields the ability to evaluate avionics sensitivities 
to existing navigation errors and hypothetical or proposed 
navigation errors. Ability to vary the error content and 
magnitude is significant during navigation system error 
budget analysis and airspace planning. 

The simulation flexibilities described above form 
the basis of the simulation test plan design. Of primary 
interest in a preliminary avionics sensitivity analysis 
are parametric values providing the greatest amount of 
useful information. During airspace planning or guidance 
algorithm development it is essential to consider the 
minimum capabilities of the navigation system. This can 
be accomplished by generating data for the worst case 
parameter values. As stated previously, one of the obj- 
ectives is to generate data for airspace planning and 
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guidance algorithm development. Hence, the desire for 
establishing worst case parameter values becomes a motivat- 
ing factor in the test plan design. 

Consideration of worst-case parameter values (flight 
profiles, geometries, error sources) in the test plan pro- 
vides data necessary to answer questions regarding air- 
space requirements, terminal area route design and RNAV/ 

MLS transition control law requirements . Airspace allot- 
ments must be made for safe separation assurance between 
aircraft equipped with minimum navigation systems suscept- 
ible to worst case error content and magnitudes . The 
subsequent route design must also provide the capability 
for the safe and expeditious passage of aircraft from 
the enroute RNAV system to final approach and landing. 

The route design must provide adequate length to execute 
maneuvers to null out any undesirable RNAV delivery errors 
detected by the more accurate MLS system. There may also 
exist some speculation as to the requirement for a guidance 
policy and/or algorithm to null out the RNAV delivery errors 
If the worst case geometries do not require implementation 
of specialized transition guidance algorithms then it is 
obvious that less than worst case geometries would not 
require additional algorithm development. 

Additionally, it is desirable to consider other more 
typical geometries to broaden the scope and, hence, the 
usefulness of the data. The additional parametric values 
of the test plan will provide greater insight into air- 
craft navigation system performance in the transition region 
and thereby aid in the route design problem and in the 
determination of the flexibility of various geometries and 
flight profile configurations. The broader data base is 
also useful for any subsequent flight test designs. 
Generally, flight test profiles are based on knowledge 
gained through analytical means or good engineering judg- 
ment. With knowledge gained from the data base generated 
by simulation techniques, more judicious selection of 
actual flight test profiles is possible. 

Implementation of the simulation is initiated upon 
completion of the test plan design. Each of the flight 
profiles is used to generate the nominal trajectory about 
which the covariance matrix is propagated. The navigation 
aid locations, because of the geometry dependence of the 
navigation errors, and the error parameters drive the error 
covariance about the nominal. Execution of the test plan 
using the ensemble simulation provides the data necessary 
for the avionics sensitivity analysis data base. 
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3.3 EXPERIMENT-SIMULATION PLAN 

The parameters to be considered for the avionics 
sensitivity study include: 

(a) Geometric parameters 

vertical profile 
horizontal profile 
VOR/DME facility location 
DME/DME facility location 
flight path/runway orientation 

(b) Navigation systems 

airborne equipment 
ground navigation aids 

(c) Error source parameters 

type (bias, correlated, random noise) 
magnitude 

The parameters to be held constant for this study are the 
vertical profile (3° CTOL approach paths) and airborne 
equipment (INS) to be compatible with the TCV aircraft 
avionics. The following test plan describes the variations 
considered for the remaining parameters. 

Two sets of horizontal profiles will be considered. 

The first set consists of the ICAO demonstration profiles 
as shown in Figures 3.1 and 3.2. Tables 3.3 and 3.4 pre- 
sent the waypoint data for each profile, respectively. 

Three hypothetical profiles are included in the second set 
in an attempt to provide more generality to the results and 
to examine worst-case situations. The three hypothetical 
profiles include a straight- in approach, a 90° -turn on to 
final and a 180°-turn on to final. 

The ground navigation aids being considered for this 
study are DME/DME and VOR/DME. The location of the ground 
facilities have been selected for the ICAO demonstration 
flights to represent those used during the demonstration. 
The demonstration flights are to be executed in the INS/ 
DME/DME (IDD) mode only. The airborne system has an 
autotune feature which selects ground facilities according 
to specific criteria. The primary facility is fixed (in 
this case, ACY) and the secondary facility is autotuned 
from Millville, Coyle or Sea Isle. For the simulation 
tests the facility providing the best geometry will be 
selected. 
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Figure 3.1 ICAO Demonstration Profile 4 
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Table 3. 3 

Altitude and GPI04 Based Upon A Calculated 
GPIP Opposite ELI Antenna 


ICAO DEMONSTRATION WAYPOINT LOCATIONS FOR PROFILES (STAR4AC043) 

WAYPOINT 

LATITUDE 

— 

LONGITUDE 

ALTITUDE, MSL 


39° 30' 31.00" 

74° 32' 42.38" 

Climbi ng 

3AC02 

39° 32' 44.62" 

74° 33' 35.09" 



39° 35' 06.59" 

74° 39' 19.50" 

4,000' 


39° 32' 00.20" 

74° 43' 46.25" 

4,000' 


39° 30' 20.70" 

74° 42' 54.6 " 

4,000' 


39° 26' 45.79" 

74° 41' 02.69" 

2,771 ' 


39° 24' 26.61" 

74° 39' 50.17" 

1 ,975' 


39° 24' 19.56" 

74° 36' 57.14" 

1,083' 

GP104 

39° 27' 07.35" 

74° 35' 02.05" 

66' 


DME 

, ARC CENTERS 


POINT 

LATITUDE 

LONGITUDE 

RADIUS 

A 

39° 3V 17.42" 

74° 34' 3^.93" 

10,000' 

B 

39° 32' 55.79" 

74° 40' 49.32" 

15,000' 

C 

39° 24' 54.36" 

74° 38' 21.55" 

7,500' 
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Table 3'. 4 


Altitude and GPI04 Based Upon A Calculated 
GPIP Opposite Ell Antenna 


ICAO DEMONSTRATION WAYPOINT LOCATIONS FOR PROFILE 3 (STAR 4AC043) 

WAYPOINT 

LATITUDE 

LONGITUDE 

ALTITUDE, MSL 

lACOl 

39° 30' 

31.00" 

74° 32' 

42.38" 

Climbing 

2AC02 

39° 32' 

14.05" 

74° 36' 

19.34" 


DDSOl 

39° 26' 

10.86" 

74° 41' 

47.36" 

4,000' 

DDSO? 

39° 25' 

23. 82" 

74° 42' 

57.17" 

4,000' 

DDS03 

39° 22' 

41 .51" 

74° 49' 

47.46" 

4,000' 

DDS04- 

39° 19' 

22.85" 

74° 45' 

24.86" 

4,000' 

DDS05 

39° T9-' 

57^@IT‘- 

74° 45’ 

01.10" 

4,000' 

Acsn 

39° 22' 

00.44" 

74° 43' 

36.82" 

3,273' 

AC522 

39° 23' 

46.48" 

74° 40' 

24.08" 

2,636' 


39° 24' 

17.12" 

74" 40' 

14.80" 

2,019" 


39P 24' 

03.02" 

74° 39' 

40.61" 

1,860' 

AC525 

39° 23' 

48.93" 

74° 39' 

06.4 " 

1 ,702' 


39° 24' 

19.56" 

74° 36' 

57.14" 

1,083' 

■■ 

39° 27' 

07.35" 

74° 35' 

02.05" 

66' 



DME ARC 

CENTERS 



POINT 

LATITUDE 

LONGITUDE 

RADIUS 

A 

39° 31' 

17.42" 

74° 34' 

34.93" 

10,000' 

B 

39° 20' 

32.37" 

74° 48' 

13.82" 

15,000' 

C 

39° 24' 

54.36" 

74° 38' 

21.55" 

7,500' 

D 

39° 27' 

35.89" 

74° 44' 

24.08" 

15,000' 

E 

39° 27' 

11.65" 

74° 40' 

59.72" 
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For the hypothetical profiles, both the DME/DME and 
VOR/DME navigation modes are assumed. The ground facility 
locations are selected to provide worst-case situations. 

The most adverse geometries are selected since they 
provide the most useful information by indicating an upper 
bound to the expected error magnitudes. For the DME/DME 
navigation modes the facilities are located such that the 
relative bearing between the two stations measured at the 
aircraft is 45° or 135° (the threshold value for the air- 
borne RNAV system) at the intersection of the flight path 
and MLS coverage boundary. When the relative bearing 
reaches these threshold values the aircraft position errors 
due to DME errors is maximum. 


The three horizontal profiles and DME facility loca- 
tions for the DME/DME mode are shown in Figures 3.3 through 
3.5. For each profile two sets of DME facility locations 
are assumed denotes by (aj^,a 2 ) and (bQ^,b 2 ) in the figures. 

The facilities (a^,a 2 ) will provide the worst-case along- 

track and the facilities (bj^,b 2 ) will provide the worst-case 
cross-track. 


The VOR/DME facility locations are shown in Figures 
3.6 through 3.8. The facility at a, in each case, rep- 
resents the worst-case along-track, due to VOR bearing 
errors. The facility at b represents the worst-case 
cross- track error due to VOR bearing errors. The worst- 
case along-track error is of interest, since for 3D RNAV 
systems, this error propagates into the vertical. All of 
the geometries are presented explicitly in Appendix B. 


The VOR and DME error sources are assumed to consist 
of random noise and bias . The magnitude of these errors 
are assumed to be consistent with those specified in Advisory 
Circular AC90-45A [10]. Although the Advisory Circular does 
not delineate the amount of bias or noise, it is assumed 
for the study that this is a variable for VOR, in the case 
of VOR/DME, and for DME in the case of DME/DME. When studying 
VOR/DME, the DME noise content does not have the significant 
effect of VOR noise content, particularly VOR bias. Hence, 
for VOR/DME the error content and magnitudes to be considered 
are : 


(a) VOR bias 
DME bias 


0.0°; VOR noise = 3.65° 

660 meters; DME noise = 660 meters; 


(b) VOR bias 
DME bias 


2.51°; VOR noise = 2.57° 

660 meters; DME noise = 660 meters; 


(c) VOR bias 
DME bias 


3.56°; VOR noise = 0.0° 

660 meters; VOR noise = 660 meters. 
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Figure 3. .6 Straight- In Approach for VOR/DME Navigation 
Mode 




Figure 3.7 90° Turn to Final VOR/DME Navigation Mode 





Figure 3.8 180° Turn to Final for VOR/DME Navigation 

Mode 
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Similarly, for the DME/DME mode the cases to be considered 
are ; 


(a) 

DME 

bias = 0.0 

meters ; 

DME 

noise = 945 

meters 

(b) 

DME 

bias = 660 

meters ; 

DME 

noise = 660 

meters 

(c) 

DME 

bias = 945 

meters ; 

DME 

noise = 0.0 

nmi 


The AC90-45A values are (for both ground and airborne) 
VOR = 3.56°; DME = 945 meters 
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IV. ANALYTICAL RESULTS 


4.1 INTRODUCTION 

As indicated previously, judicious selection of the 
experimental test parameters (VOR/DME location, DME/DME 
location, flight path geometry, runway orientation [MLS 
antenna orientation] and error source magnitudes) will 
provide the capability to generate an extensive data base. 

For each of these configurations, the avionics sensitivity 
to navigation errors for the 737-TCV aircraft were evaluated. 
The data base for this study consists of 36 cases, 18 DME/ 

DME transitions to MLS and 18 VOR/DME transitions to MLS. 

Meaningful presentation of the results becomes a dif- 
ficult task in light of the scope of the data base. The 
primary quantities of interest are the cross-track and 
vertical deviation time histories through the RNAV/MLS 
transition as a function of the various study parameters. 
These time histories are presented in the form of plots. 

An attempt has also been made to tabulate several of the 
more important quantities, such as, the RNAV deviations at 
the MLS boundary as a function of the study parameters. 

The results section outline is as follows. The simul- 
ation test conditions are presented in Section 4.2. These 
conditions essentially establish the ground rules for the 
subsequent experimentation program. Section 4.3 presents an 
explanation of the contents of each of the plots and on how 
each plot should be read. This explanation is performed 
via an example. The avionics sensitivity analysis results 
are presented in Section 4.4. This section contains the 
results of the entire sensitivity study. Demonstration 
results are shown in Section 4.5 which contains results 
for the ICAO demonstration flights. A brief summary of 
results is presented in Section 4.6 for quick reference. 

4.2 SIMULATION TEST CONDITIONS 

The purpose of this subsection is to identify those 
conditions established for the simulation experiments. 

These conditions relate primarily to the initialization of 
the covariance matrix. Specifically, the initialization 
of the covariance elements relating to the INS, altimeter 
and position estimates. Reinitialization of the covariance 
matrix at MLS transition is also discussed where it is as- 
sumed that the avionics switches directly from the RNAV mode 
to the MLS mode. 
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The INS initial position error estimates are dependent 
on the most recent radio navigation update. For example, 
an INS- equipped aircraft traversing the ocean without an 
update will experience a drift error of approximately 7413 
meters assuming a drift rate of 1853 meters/hr and a flight 
time of 4 hours. This drift error is reduced, through 
external position fixes, to an amount equivalent to the 
measurement errors associated with the external source. 

For this study it is assumed, depending on the mode of 
operation, that the INS has been updated using VOR/DME or 
DME/DME. It is further assumed that the INS initial posi- 
tion errors are equivalent to the VOR/DME or DME/DME 
position errors at the initial position of the simulation 
test runs. Hypothetically, the INS could have been updated 
at more favorable geometries thereby decreasing the initial 
INS error estimate. However, for this study the motivation 
is to determine upper bounds, hence, it is more desirable 
to examine worst-case conditions. The magnitude of the 
initial INS errors is therefore determined using AC 90-45A 
[10] techniques to establish VOR and DME derived position 
errors for a given aircraft location relative to the facil- 
ity . 


For example, consider the initial geometry shown in 
Figure 4.1. An aircraft with a range of 20.7 km from the 
facility on a bearing of 154° will have VOR/DME- induced 
position errors, for an aircraft heading of 0°, of 730 
meters along-track and 1170 meters cross-track [exclusive 
of flight technical and computer errors) . Since the air- 
craft course for this example is due north the along-track 
error becomes the initial INS north position error and the 
cross-track error becomes the initial INS east position 
error. 

The other elements of the INS are represented by the 
following quantities . The initial INS velocity error is 
considered to be 0.05 meters/sec. The gyro drifts, con- 
sidered as biases because of the short flight times, are 
0.01 deg/hr, the accelerometer errors, also assumed biases 
for this study, are 10"4g, where g is the gravitational 
constant. The bias assumptions are valid in light of the length 
of the correlation times of 5 hours for the gyro drift 
and 10 hours for the accelerometer relative to the flight 
times on the order of 5 to 10 minutes. The platform mis- 
alignment for this study is taken to be 0.005°. 

The altimeter model for this study is assumed to con- 
sist of a bias of 30.5 meters and random noise of 7.6 meters. 
The altitude measurement is passed through a first order 
filter with a time constant of 2 seconds to simulate alti- 
meter lag. The initial covariance for the altimeter is 
taken to be 30.5 meters, consistent with the altimeter bias. 
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The initial aircraft cross-track error is arbitrarily assumed 
to be 927 meters. 

The covariance matrix is propagated through the RNAV 
region using the initial conditions specified above. At 
the MLS boundary the covariance is reinitialized to rep- 
resent the fact that the RNAV measurements and filters are 
replaced by the appropriate MLS measurements and filters. 

The elements of the covariance matrix common to both the 
RNAV and MLS simulations remain unaltered through transition. 
These include the aircraft dynamics, autopilots and vertical 
complementary filter. The diagonal terms of the MLS comple- 
mentary filter are set at the diagonal values of the RNAV 
complementary filter. The MLS accelerometer and platform 
alignment inputs (from the AHRU) are assumed to be the same 
as the INS. The MLS measurement biases are set in accordance 
with Table 3.1. 

To conserve computer time in order to examine as many 
test cases as possible, the length of the simulation tests 
was shortened. More specifically, the simulation was termin- 
ated after the covariance reached a steady state (except for 
geometry variations) after transition into the MLS. This 
still provided the most useful information for the current 
study on RNAV/MLS transition behavior. 

4.3 PRESENTATION OF RESULTS 

The results of the sensitivity analysis will be displayed 
using 20 plots. Ten plots will contain the cross-track devi- 
ation time history and the other 10 plots will contain the 
vertical deviation time histories. Each corresponding cross- 
track and vertical deviation time history plot is associated 
with a particular geometry. The geometric configurations 
consist of three flight profiles (straight-in approach, 

90“ turn- to- final , and 180® turn- to- final) and two ground 
facility locations for both VOR/DME and DME/DME. 

Figure 4.2 shows an example cross-track deviation time 
history plot. Each of the other plots presented generally 
have the same information content. Three curves are presented 
on each plot depicting the avionics sensitivity to navigation 
sensor errors. For this particular example the curves cor- 
respond to : (1) a VOR bias of 3.56° and no noise component 

(bias), (2) a VOR bias and noise component of 2.51° (bias/ 
noise), and (3) a VOR noise component of 3.56° without bias 
(noise). The DME errors for all three cases is DME bias and 
noise components of 660 meters. On the time axis are indi- 
cated the two points defining the straight and level portion 
to allow initial condition transients to dissipate (between 
waypoints and WP 2 ) and the point of MLS transition. 
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Figure 4.2 Example Result Plot for Illustration 



Hence, each plot will reveal sensitivities to error para- 
meter variations. 

Geometric sensitivities can be determined by comparing 
various plots. To facilitate this comparison, Table 4.1 
has been established as a cross-reference guide. The "figure 
numbers" in column one correspond to a particular figure in 
the results. The geometry and its corresponding figure 
number are shown in the third and fourth columns. The final 
two columns indicate the mode of operation, VOR/DME or DME/ 
DME, and the facilities being used to provide navigation 
support (the letters correspond to the letters in the ap- 
propriate geometric figure) . 

4.4 SENSITIVITY ANALYSIS 

The presentation of the sensitivity results is broken 
down into two subsections for clarity. The first subsection 
presents the VOR/DME results and the second the DME/DME re- 
sults . The motivation stems from the fact that most of the 
DME/DME cases can be represented by a VOR/DME case. For ex- 
ample, the DME/DME delivery errors, for facilities and B 2 

(worst-case cross-track) and DME bias and noise of 465 meters, 
are 355 meters in cross -track and 38 meters in vertical devi- 
ation. The comparable VOR/DME case is facility A (worst- 
case along-track) with a VOR bias of 0° and noise of 3.56°. 

The RNAV delivery errors associated with this case are 394 
meters in cross-track and 33 meters in vertical deviation. 

To conserve computer time it is assumed that these two cases 
would have the same transient response in the MLS region 
since they represent almost identical initial conditions at 
the MLS boundary. Table 4.2 indicates the RNAV delivery 
errors for each geometry tested and indicates which of the 
DME/DME case was selected to be represented by which VOR/DME 
case. For example, the cross-track transient response of 
DME/DME case h is represented by the cross-track transient 
response of VOR/DME case g. 

4.4.1 VOR/DME 


The cross-track and vertical deviation time histories 
for the VOR/DME RNAV test cases are shown in Figures 4.3 
through 4.14. Figures 4.3 and 4.4 show the straight-in 
approach for VOR/DME facility A (worst-case along-track). 
During the RNAV mode the observation can be made that the 
initial condition transients have decayed during the straight 
and level portion of the flight path. For the vertical devi- 
ation pitch over is easily recognized by the vertical channel 
response to the along-track feeding in (Figure 4.4) which is 
typical of 3D RNAV systems. Also of special interest is the 
sensitivity of the RNAV filter to sensor error content. 
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Table 4.1 

Sensitivity Analysis Results Plots Cross-Reference Table 


RESULTS 

FIGURE 

NUMBER 

TIME 

HISTORY 

GEOMETRY 

GEOMETRY 

FIGURE 

NUMBER 

NAVIGATION 

MODE 

FACILITY(IES)* 

4.3 

CTD 

Straight-In 

3.6 

VOR/DME 

A 

4.4 

VD 

Straight- In 

3.6 

VOR/DME 

A 

4.5 

CTD 

90° Turn 

3.7 

VOR/DME 

A 

4.6 

VD 

90° Turn 

3.7 

VOR/DME 

A 

4.7 

CTD 

180° Turn 

3.8 

VOR/DME 

A 

4.8 

VD 

180° Turn 

3.8 

VOR/DME 

A 

4.9 


Straight- In 

3.6 

VOR/DME 

B 

4.10 


Straight-In 

3.6 

VOR/DME 

B 

4.11 


90° Turn 

3.7 

VOR/DME 

B 

4.12 


90° Turn 

3.7 

VOR/DME 

B 

4.13 

CTD 

180° Turn 

3.8 

VOR/DME 

B 

4.14 

VD 

180° Turn 

3.8 

VOR/DME 

B 

4.15 

CTD 

Straight- In 

3.3 

DME/DME 

A^ , A2 

4.16 

VD 

Straight-In 

3.3 

DME/DME 

Ap A2 

4.17 

CTD 

90° Turn 

3.4 

DME/DME 

A^, 

4.18 

VD 

90° Turn 

3.4 

DME/DME 

A^. A 2 

4.19 

CTD 

180° Turn 

3.5 

DME/DME 

A-| ) A2 

4.20 

VD 

180° Turn 

3.5 

DME/DME 

A.|, A 2 

4.21 

CTD 

— 

3.3 

DME/DME 

Bi , B2 

4.22 

VD 

— 

3.3 

DME/DME 

Bi* B2 


* Facility A; 
Facility B: 


Worst-Case Along-Track 
Worst-Case Cross-Track 
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Table 4.2 


MLS Transition Cases Executed to Complete Simulation Test Plan 


CASE 

MODE 

FACILITY 

ERROR 

RNAV' DELIVE 

RY ERROR 

C0RRESP0^ 

DING CASE 

CROSS-TRACK 

(METERS) 

VERTICAL 

(METERS) 

■IMiM 

VERTICAL 

a 

VOR/DME 

A 

Bias 

407 

64.3 

* 

★ 

b 


(Along-Track) 

Bias/Noise 

400 

50.9 

★ 

* 

c 



Noise 

394 

32.6 

* 

* 

d 


B 

Bias 

1057 

37.2 

★ 

★ 

e 


(Cross-Track) 

Bias/Noise 

755 

37.2 

* 

★ 

f 



Noise 

169 

37.2 

★ 

★ 

g 

DME/DME 

Ai 

Bias 

323 

46.3 

* 

* 

h 


(Along-Track) 

Bias/Noise 

245 

39.6 

g 

d 

i 



Noise 

126 

31.7 

f 

c 

j 


B] .Bg 

Bias 

486 

43.0 

a 

d 

k 


(Cross-Track) 

Bias/Noise 

355 

37.8 

c 

d 

1 



Noise 

131 

31.7 

f 

c 


* MLS transitions were executed for these cases. 
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Figure 4.3 Cross-Track Deviation for Straight- In, VOR/DME , Facility A 



VCRTICflL OEVIflTIdN (hETCRS) 








CRBSS-TRRCK DEVinTlBN tHETERS) 



0 SO 100 150 200 250 300 35 


Ttne (SECBNDS) 

Figure 4.5 Cross-Track Deviation for 90° Turn to Final, VOR/DME, Facility A 
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Figure 4.6 Vertical Deviation for 90® Turn to Final, VOR/DME, Facility A 






CR8SS-TRRCK OEVlRTieN (nETERS) 



0 50 100 150 200 250 300 350 

Tins (SEC8N0S) 


Figure 4.7 Cross-Track Deviation for 180° Turn to Final, VOR/DME, Facility A 
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Figure 4.8 Vertical Deviation for 180® Turn to Final, VOR/DME, Facility A 
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Figure 4.11 Cross-Track Deviation for 90® Turn to Final, VOR/DME, Facility B 
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Figure 4.12 Vertical Deviation for 90® Turn to Final, VOR/DME, Facility B 
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Figure 4.13 Cross-Track Deviation for 180® Turn to Final, VOR/DME, Facility B 
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Figure 4.14 Vertical Deviation for 180® Turn to Final, VOR/DME, Facility B 



With pure random noise (high frequency) the filter displays 
excelelnt performance. Addition of bias into the signal 
content significantly degrades filter performance. 

During MLS transition it can be observed that the initial 
transients are significant. This is particularly true of the 
vertical channel (note also that the vertical channel has a 
significant transient at pitch over initiation). Several 
offline studies indicated that the magnitude of the response 
was sensitive to filter initialization at the onset of MLS 
transition. However, initialization did not impact the 
settling time. 

At MLS transition the filters are subjected to a step 
input (bias) which can be significant. The filter and auto- 
pilot coupling during transition was found to be the primary 
cause for the subsequent transient responses to this input bias 
One factor not included in the ensemble simulation because of 
the linear characteristics of this technique, is the control 
bound normally imposed on a control system. This will reduce 
the magnitude of the transient and will probably lengthen the 
settling time. However, this needs to be verified through 
nonlinear simulation techniques or through describing func- 
tion techniques in the ensemble simulation. 

Figures 4.5 and 4.6 show the cross-track and vertical 
deviation time histories respectively for the 90° turn-to- 
final approach. Similarly, Figures 4.7 and 4.8 show the 
results for the 180° turn- to- final approach. Note that 
the variation of the magnitude of the errors and the set- 
tling times of the position errors in the MLS region vary 
as a function of geometry. This is primarily due to the 
range dependence of the azimuth and elevation- induced posi- 
tion errors. Hence, the MLS geometry effects impacts both 
the magnitude and the settling times. This implies that the 
route design and MLS coverage region must provide sufficient 
airspace to null the RNAV delivery errors such that air- 
craft stabilization is achieved prior to final approach and 
landing . 

Figures 4.9 through 4.14 display the cross-track and 
vertical deviation time histories for the VOR/DME RNAV cases 
with the worst-case cross-track delivery errors. Note that 
the vertical time histories are the same for all error para- 
meter values. This is due to the fact that the along-track 
error is driven by the DMF errors only which, for the VOR/ 

DME cases, are not varied. Similar transient behavior is 
noted with these cases as with the worst-case along-track 
error cases discussed previously. 

An additional peculiarity of interest is the transient 
response due to increased delivery errors as observed in 
Figures 4.9, 4.11 and 4.13. For larger cross-track errors 
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the MLS transients display additional oscillations. Par- 
ticularly in the 180° turn- to- final (Figure 4.13), several 
additional oscillations are observed (these are probably 
geometry- induced since the aircraft/runway orientation is 
continually changing throughout the turn) . For these cases 
the settling time is assumed to occur after the second 
oscillation since the response appears to be approaching 
steady state at this point. 

In Figures 4.10, 4.12, and 4.14 the vertical channel 
transients are observed to be significant. Again it is 
partially attributable to initial condition transients and 
to f ilter/autopilot coupling. 

4.4.2 DME/DME 

Figures 4.15 through 4.22 indicate the results for the 
DME/DME updated INS RNAV simulation. Only Figures 4.15 
through 4.20 show transition into MLS for the reasons stated 
previously. Generally, the DME/DME- induced RNAV delivery 
errors are less than the VOR/DME- induced delivery errors. 

For this reason the transients tend to be more well behaved 
during transition. The improved navigation performance 
characteristics of the DME/DME system (relative to the VOR/ 
DME system) is also observed in the RNAV portion. 

Figures 4.21 and 4.22 indicate the worst-case cross- 
track tests for DME/DME. To determine the MLS transient 
response the corresponding VOR/DME transient is consulted 
as indicated in Table 4.2. 

4.5 ICAO DEMONSTRATION FLIGHTS 

For illustrative purposes the 130° turn- to- final and 
"S-turn" approach utilized during the ICAO demonstration 
flights were simulated. The DME/DME facilities selected 
for each case were Sea Isle and Millville. These were 
chosen since it was difficult to determine a facility pair, 
using the Atlantic City VORTAC as primary, that provided 
a favorable geometry. The error associated with these 
facilities was assumed to be a DME bias of 350 meters and 
DME noise of 650 meters. The waypoints associated with 
the ICAO demonstration flight profiles are indicated on the 
plots shown in Figures 4.23 through 4.26. These correspond 
to the profiles shown in Figures 3.1 and 3.2. 

4.6 SUMMARY OF RESULTS 

The significant results of the sensitivity analysis 
are presented here as a quick-look summary. The summary 
is given in Table 4.3. The table includes the RNAV delivery 
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Figure 4.15 Cross-Track Deviation for Straight-In, DME/DME, Facilities A 
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Figure 4.17 Cross-Track Deviation for 90° Turn to Final, DME/DME , Fa-cilities A 
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Figure 4.18 Vertical Deviation for 90“ Turn to Final, DME/DME, Facilities A 
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Figure 4.20 Vertical Deviation for 180® Turn to Final, DME/DME, Facilities A 
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Figure 4.24 Vertical Deviation Time History for S-Turn ICAO Profile 




CROSS TRUCK DEV1RT10N (METERS) 



Figure 4.25 Cross-Track Deviation Time History for 130° Turn ICAO Profile 
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Figure 4.26 Vertical Deviation Time History for 130° Turn ICAO Profile 



Table 4.3 
Summary Results 


MODE 

FACILITY 

ERROR 

RNAV DELIVERY 

MLS SETTLING TIMES (MINUTES) 

ERR0RS(M 

ETERS) 

STRAIGHT-IN 

90° TURN 

180° TURN 1 

CROSS- 

TRACK 

VERTI- 

CAL 

CROSS- 

TRACK 

VERTI- 

CAL 

CROSS- 

TRACK 

VERTI- 

CAL 

CROSS- 

TRACK 

VERTI- 

CAL 

VOR/DME 

Worst-Case 

Bias 

407 

64.3' 

0.95 

2.33 

0.95 

1.63 

0.95 

1.54 


Along-Track 

Bias/Noise 

400 

50.9 

0.95 

2.33 

0.95 

1.63 

0.95 

1.54 



Noise 

394 

32.6 

0.95 

2.33 

0.95 

1.63 

0.95 

0.83 


Worst-Case 

Bias 

1057 

37.2 

1.77 

0.95 

1.48 

0.95 

* 

0.85 


Cross-Track 

Bias/Noise 

755 

37.2 

1.67 

0.95 

1.35 

0.95 

if 

0.85 

1 


Noise 

169 

37.2 

0.95 

0.95 

0.95 

0.95 

0.90 

0.85 

DME/DME 

Worst-Case 

Bias 

323 

46.3 

0.95 

* 

0.95 

0.95 

0.95 

0.95 


Along-Track 

Bias/Noise 

245 

39.6 

0.95 

0.95 1 

0.95 

0.95 

0.95 

0.95 



Noise 

126 

31.7 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 


Worst-Case 

Bias 

486 

43.0 

0.95 

* 

0.95 

0.95 

0.95 

0.95 


Cross-Track 

Bias/Noise 

355 

37.8 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 


1 

Noi se 

131 

1 

31.7 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 


* Oscillations remain at end of simulation run (2.5 minutes after MLS transition) 




errors at' the MLS boundary for each geometry and error magni- 
tude considered in the study. The associated settling time 
was estimated from the plots for each of the MLS geometries 
(flight profile/runway configuration) . 

The settling times reflect the impact of two contribut- 
ing factors. The first is the magnitude of the RNA"V 
delivery errors and the second is the magnitude of the 
range-dependent MLS errors (azimuth and elevation) . The 
latter factor implies the need for range-dependent filter 
gains to minimize this effect. 
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V. CONCLUSIONS 


The primary objective of this study was to generate 
a data base of position errors for the TCV aircraft in the 
RNAV/MLS transition region. From the data base an analysis 
is performed to evaluate TCV avionics sensitivity to flight 
profile, sensor errors, ground facility location and runway 
orientation (MLS antenna orientation) . The avionics sensi- 
tivity analysis provides insight into such issues as RNAV/ 

MLS transition feasibility (for various geometries and flight 
profiles), transition guidance requirements and MLS airspace 
and operational procedure requirements. The data base is 
also useful for subsequent flight test planning. 

Generally, the data base consists of cross-track and 
vertical deviation time histories from the RNAV region through 
MLS transition. The RNAV system performance is best assessed 
by examination of the delivery errors at the MLS boundary. 

For the DME/DME updated inertial RNAV system the cross-track 
delivery error varies from 126 meters to 486 meters. The 
vertical delivery error ranges from 31.7 meters to 46.3 
meters. For the VOR/DME updated system the cross- track 
delivery error ranges from 169 meters to 1057 meters and the 
vertical delivery error ranges from 32.6 meters to 64.3 meters. 
The improved navigation system performance in the DME/DME mode 
as opposed to the VOR/DME mode is obvious from these ranges. 

The error source magnitudes for this study were assumed 
to be consistent with the quantities specified in AC 90-45A 
[10]. This document, however, does not specify the relative 
content of random bias and noise value for each navigation 
source. The amount of bias and noise was varied parametric- 
ally to assess avionics sensitivity to signal noise content. 

It was found that the complementary filter performance is 
sensitive to signal content in that delivery error magnitude 
increased significantly with bias content increase. For 
example, for the VOR/DME RNAV system the cross- track delivery 
error increased from 169 meters for all noise and no bias to 
1057 meters for all bias and no noise. Similarly, the vertical 
delivery error ranges from 32.6 meters for all noise and no 
bias to 64.3 meters for the all-bias case. The corresponding 
DME/DME delivery increases were from 131 meters to 486 meters 
in cross-track error and from 31.7 meters to 46.3 meters in 
vertical error. 

An area of concern in RNAV/MLS transition is the require- 
ment for additional guidance during transition. This require- 
ment stems from the need to null out RNAV delivery errors upon 
transition to MLS in an expeditious and smooth manner. The 


75 



requirement will depend on the magnitude of the delivery 
errors. For the DME/DME updated inertial system the de- 
livery errors were considered sufficiently small not to 
warrant additional guidance transition. The magnitude 
of the delivery errors for the VOR/DME updated inertial 
system were considered to be sufficiently significant to 
require additional transition guidance^ 

The MLS complementary filters were observed to dis- 
play significant transients immediately upon transition. 

This behavior was attributed to two factors. First of 
all, an extensive amount of cross - coupling exists between 
the filters and the autopilot/aircraf t dynamics. Secondly, 
the magnitude of the transients was found to be sensitive 
to filter initialization. The transient settling time was 
found to be primarily a function of the system time con- 
stants and input magnitudes (RNAV delivery errors and range 
dependent MLS errors). To minimize initialization transients 
the filters were initialized with the RNAV filter estimates 
at transition. The transient response magnitude was thereby 
reduced. 

As indicated, however, the magnitude of the RNAV 
delivery errors and the magnitude of the MLS sensor errors 
(especially the range- dependent azimuth and elevation errors) 
have a significant impact on the settling time. For certain 
flight profiles and delivery errors, additional oscillations 
are induced thereby increasing the settling time. The set- 
tling times for the cases examined during this study varied 
from 0.85 minutes to 2.33 minutes (although several cases had 
not settled at the end of the simulation run) . 

These settling times have a significant impact on air- 
space requirements and operational procedures in the MLS 
environment. This is particularly true for approach speeds 
of 260 km/hr typical of TCV~737 type aircraft. Hence, the 
route design and the alloted airspace must be adequate to 
provide sufficient space for aircraft stabilization prior 
to final approach and landing. 

The feasibility of a particular flight profile and 
runway/flight path configuration to support RNAV/MLS transi- 
tion is dependent on the transient response associated with 
that particular configuration. The transient response has 
potential improvements available through appropriate filter 
initialization. Other modifications are possible filter 
design changes effective immediately after transition. This 
may include some sort of blending of navigation signals or 
implementing range- dependent gains. 
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VI . RECOMMENDATIONS 


While meeting the objectives of this study, several issues 
were uncovered which require further analysis and evaluation. 
The issues stem primarily from the behavior of the transient 
response during MLS transition. The current study addressed 
rather highly sophisticated RNAV systems, INS updated with 
DME/DME or VOR/DME. It was observed that the MLS transition 
transient was dependent on RNAV delivery errors. Of interest 
would be to evaluate less sophisticated RNAV systems which 
would yield larger delivery errors. Such systems would 
include air data, updated with VOR/DME and DME/DME and course- 
line RNAV computers . 

The current study examined only a limited data base scope. 
The scope should be expanded to include additional geometries 
and flight profiles and additional error source variations 
such as MLS sensor errors. The broadened data base would 
provide greater insight into the evaluation of MLS airspace 
requirements and operational procedures. This, in turn, 
impacts the route design and the requirement for expanded MLS 
coverage (for example, an assessment of whether MLS azimuth 
coverage should be 20°, 40° or 60°). 

The transient behavior has indicated additional research 
in the area of MLS filter design. The dependence of the tran- 
sient (magnitude and settling time) on the magnitude of de- 
livery errors and subsequent flight path/runway geometry 
indicates a potential need for time or range dependent gain 
filter or filter blending techniques at MLS transition. 
Evaluation of potential filter modifications at transition 
can be accomplished through the simulation techniques des- 
cribed in this report. 

Of current interest in the navigation community is the 
evaluation of proposed navigation systems such as LORAN-C, 
Omega, and GPS. These are future systems in the same sense 
as the MLS. It would be of interest to evaluate the tran- 
sition from these systems to MLS in a similar manner as 
performed during this study. 

The final recommendation is that maximum use be made of 
the covariance propagation simulation for subsequent flight 
test design. It is possible, through implementation of this 
technique, to evaluate candidate flight test profiles and 
identify additional profiles of interest in a cost effective 
manner. 
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APPENDIX A 


RNAV/MLS ENSEMBLE SIMULATION 


A.l INTRODUCTION 

The purpose of this appendix is to present the analy- 
tical models of the RNAV and MLS simulations utilized in 
the avionics sensitivity analysis. The models are de- 
rived from information gathered from the NASA Langley 
Research Center. Unfortunately, the necessary inform- 
ation did not exist in a concise form and it became 
necessary to compile the required analytical expressions 
from various documents and computer listings. Since it 
is difficult to ascertain which documents the various models 
were extracted from it is impossible to furnish a list of 
references. Hopefully, the data gathered in this appendix 
will be useful as a future reference source. 

For the sake of convenience the models are presented 
in conjunction with the block diagrams shown in Figures 
2.1 and 2.2. These figures are repeated in the appendix 
for easy reference. 

A. 2 RNAV SIMULATION 

The VOR/DME RNAV simulation block diagram is shown 
in Figure A.l. Due to the input and output requirements 
of the various models making up the entire TCV system, 
several coordinate transformations will be presented to 
provide completeness to this discussion. 

The basic measurements of tlia TCV RNAV system are 
the range, p , bearing, 0 , the altitude, h^^ , and 

the inertial navigation system (INS) horizontal velocity 
components, Vj^ and V^ , and vertical acceleration, a^. 

The VOR/DME measurements are assumed to consist of a random 
noise and bias component (this assumption is for the cur- 
rent study only, correlated noise is readily included re- 
quiring, however, an additional state). Hence, these 
measurements are expressed as 

Pm " Pa "^p ’ 


79 



00 

o 



Figure A.l VOR/DME RNAV Simulation Block Diagram 




(A. 2) 


m 


= 0 + 
a 


+ n. 


where 



0 . 


These expressions must be linearized for the ensemble simul- 
ation. In the above expressions p and 9 are the actual 

values and bias components and (rip> ^Iq) 

the noise components. Linearizing about a nominal 

the measured and actual values become 



Pn 

+ 

^Prr. » 

m 

(A. 3) 

% = 


+ 

» 

m ’ 

(A. 4) 

Pa = 

Pn 

+ 

Ap , and 

(A. 5) 

9 = 

3 

9n 

+ 

A0 . 

(A. 6) 


Substitution of Equations (A. 3) through (A. 5) into Equations 
(A.l) and (A. 2) yields 

Apm " t>p CA. 7) 

APj^ = A0 + bg + ng (A. 8) 

Ap and A9 , the deviations of the actual values from 
the nominal (due to aircraft displacement from the nominal) 
are based on the cross-track deviation, CTD, of the aircraft. 
The along-track is not included since for 31) RNAV systems 
only an indicated and not an actual along-track deviation 
is defined. The coordinate transformation from CTD to 
Ap and A9 is given as 


Ap = CTD sin (i^-0j^), and 


(A. 9) 


A0 


CTD 

Pn 


cos 


(A. 10) 


where pj^ and 0^^ are the nominal measurement values and 
is the nominal aircraft heading. 

The altimeter model is assumed to have noise and bias 
corrupting the actual altitude. This corrupted altitude is 
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then passed through a first order filter to simulate 
altimeter delay in recording the altitude. The analy- 
tical expressions defining the altimeter model are 


where 


h^ = ~ h -t- ^ ( z b, n, ) , and 
m Th m 'h n n 


(A. 11) 
(A. 12) 


- altimeter lag, 

h^ - measured altitude, 
m 

z - actual altitude, 
bj^ - altimeter bias, and 
nj^ - altimeter noise 

Because of the large number of states required to 
fully represent the TCV RNAV system, the INS models are 
kept as simplistic as potentially possible. For this 
presentation only the final lineari 2 ed form used in the 
simulation is given. Normally, the platform drift errors 
and accelerometer errors are assumed to be correlated with 
a correlation time of 5 hours and 10 hours respectively. 


(A. 13) 
(A. 14) 


Since 

the 

flight 

times 

much 

less 

these 

errors 

model 

is 

therefore 


= - 


* 

sin 

= Pz 




= -CO 

2 5 V 

s °^N 

- (2Q, 


= CO 

s E 

+ (2Q, 

♦x 


lb - 00 li; * 

^y y 

♦y 

= -CO 

ip + 

z ^x 

00 Ip 
X ^z 


Pz) ^ ^z 


a 

z 


= 0) \b - CO li + e 
y X ^y 2 

• • • • • 

= £ = £ = a = a = ct = 0 

y z X y 2 

(time constants = 5 hrs § 10 hrs) 


a^ CA.15) 

a ip (A. 16) 
y ^z 

(A. 17) 
(A. 18) 
CA.19) 
(A. 20) 
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where Ap and A9 are outputs o£ the second order filter. 
The equations representing the second order filter are 


/N 

^N 

“ - •'i 

; X < 

+ 

*2N * 

* ’'l ’'mN ■ 

(A. 25) 

^2N 

• 

' - >'2 

^N 

+ 

*^2 ^mN . 


(A. 26) 

/N 

^E 

■ - 

^E 

+ 

’'2E *^E 

* H ’'mE . Pnd 

(A. 27) 

^2E 

= - K2 

^E 


'^2 • 


(A. 28) 


where 


*N’^E 

^2N’^2E 

6Vj^.6V^ 




filtered north and east position errors, 
filter states (north and east) , 

INS velocities (north and east) , and 
measured north and east positions. 


Integration of Equations (A. 25) and (A. 27) and an ap- 
propriate coordinate transformation provide the computed 
range, Ap , and bearing, A6 , as 

✓s ^ 

APc = cos0j^ + Xg sin6j^, and (A. 29) 


N F 

A6^ = — sine., + — cose 

C Pn n Pn 


(A. 30) 


The inputs to the vertical complementary filter are 

the measured altitude h and the vertical acceleration, 

m ’ 

a_ , from the INS. The analytical equations of the second 
order vertical complementary filter are given by 



ot. 


^l^h 


m 


^i^h 


(z.bh-rih) 



(A. 31) 
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and 


= - i- h + i- (z + b, + n, ) - z, 
me ^ 1 

where 

z^ - commanded altitude rate, 

- vertical acceleration from INS, 
bj^ - altitude bias (altimeter), 

- measured altitude (baro), 
nj^ - altimeter noise, 

- integral feedback term 

- filter time constants, and 
K*,K^ - filter gains. 


The outputs of the vertical complementary filter 
are fed directly into the vertical autopilot. The hori- 
zontal filter outputs, however, are converted into an 
estimated cross track deviation, CTD, and estimated air- 
craft heading, \jj , as 

CTD = ■ sini|j + Xg cosijj , and (A. 33) 




4 ^ = 


4 sinijj 
V 


+ 


c COSlp 

V 


(A. 34) 


where Xg, Xjj^ and Xg are filter outputs, ip is the 

nominal aircraft heading and V is the nominal aircraft 
speed. The quantities given by these two equations are 
input into the lateral autopilot according to 



(A. 35) 
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and 



where 





^y^P- 


X 


^P^g.S-i Z (A. 36) 
57.3gXp 


<(> 



X 

As. 

'P 


roll angle (input to dynamics), 
roll rate , 

time constant , 

autopilot gains , 

ground speed, 

computed cross-track error , and 
computed heading error- 


For 3D RNAV systems the along- track deviation, ATD, 
propagates into the vertical channel. As indicated in 
Figure 2.3 in the main text, an indicated along- track 
deviation will shift the vertical profile so as to pro- 
duce a vertical deviation. The indicated along-track is 
computed from the measured and filter computed range and 
bearing according to 

ATD = (Apjjj-Ap^)cos (0j^-i|j) + (A0jjj-Ae^) pj^sinCBj^-iij) (A. 37) 


where Ap^ and A0^^ are given by Equations (A. 7) and 

(A. 8), Ap^ and A0^ are given by Equations (A. 29) and 

(A. 30), Pj^ and 0j^ are the nominal range and bearing, 

and \p is the nominal aircraft heading. With this along- 
track deviation and the outputs of the second order vertical 
complementary filter, the vertical autopilot equations are 

h^ = (A. 38) 


85 



and 


^ 1 h + 1 

^2 ' T ^2 T 




09 (h^-hp 


(A. 39) 


where z is obtained from Equation (A. 39), is the de- 
sired altitude, and x is the autopilot time constant. 

The aircraft dynamics are expressed as second order 
in the lateral degree of freedom and second order in the 
vertical degree of freedom. The four equations, therefore, 
expressing the dynamics are 


CTD = V, 


(A. 40) 


= f 4>, 


(A. 41) 


z = h, , and 


>'l - ^2 


(A. 42) 
(A. 43) 


where (p is given by Equation (A. 35) and h 2 is given by 
Equation (A. 39). The remaining variables are 

CTD - cross track deviation, 

V - nominal aircraft speed, 
ip - nominal aircraft heading 
g - gravitational constant, and 
z - altitude. 

The DME/DME RNAV simulation block diagram is given in 
Figure A. 2. The primary difference between this simulation 
and the VOR/DME simulation depicted in Figure A.l are several 
coordinate transformations between various analytical models. 
The coordinate transformations are different for DME/DME 
since these measurements are in a different frame of refer- 
ence. 
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Figure A. 2 DME/DME RNAV Simulation Block Diagram 






First o£ all, the transformation from the actual 
aircraft CTD to the measurement errors requires a trans- 
formation to north and east error components. Then error 
components, along with a vertical error, are transformed 
into two DME errors . The first transformation is given by 


AXj^ = - CTD sinip, and (A. 443 

AXg = CTD cosip (A. 45) 

where 


AXj^ - north position error, 

AXg - east position error, 

CTD - actual cross track deviation, and 
- nominal aircraft heading 

The transformation from north and east errors to the measure- 
ment errors, Ap , and Apo, is given by 

1 ^ 


^1 


(A. 46) 


Ap = — Ax., + — ^ ■ AXp + 
2 p„ N p,, E 


(hj-hp 


Ah, 


CA.47) 


where 


p, - nominal slant range from reference facility to 
aircraft, 

x^ - nominal aircraft position north of reference 
facility, 

Xp - nominal aircraft position east of reference 
facility , 
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- displacement o£ secondary facility north of 
reference facility, 

y. 2 v ~ displacement of secondary facility east of 
reference facility, 

- nominal aircraft altitude with respect to 
reference facility, 

h« - altitude of secondary facility with respect 
to reference facility, and 

P 2 - nominal slant range from secondary facility 
to aircraft given by 

1/2 

^2 " (h]^-h2)^} (A. 48) 

These quantities are shown schematically in Figure A. 3. 

The two DME measurements are given as 


= Apj^ ^p ^p’ (A. 49) 

AP2^^ = AP2 + bp + np CA.50) 

where 

b - DME bias 
P 

n - DME noise 
P 

Ap^ - given in Equation (A. 46) 

Ap 2 - given in Equation (A. 47), and 

Apj^™, Ap 2 ^ - DME/DME measurement errors 

The differential equation governing bp is 

tp = 0 (A. 51) 

Two coordinate transformations remain the first trans- 
forming the DME/DME measurements into north and east comp- 
onents to be compatible with the complementary filter inputs 
and the second deriving the indicated along-track deviation 
from the north and east error components. The first 
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transformation is 



M (- 

. -a 

(A. 52) 


M ■ ( 

^N’^2N\ a - ^ . 

(A. 53^ 

where 

^ ^E ^2N 

^2E 

CA. 54) 


- PjP2 

• 

It is of interest to note that the quantity D approaches zero 
as the angle between the slant ranges approach 0° or 180°, a 
peculiarity of the DME/DME mode of navigation. In these ex- 
pressions Ap^ and Ap^ are given by 


A - A in 

Apj^ = Ap^ - 

Apj^^ , and 

(A. 55) 


Apz = Apz - 

A c 

Apz 

(A. 56) 

where 

A ni . m 
Ap^ , Apz 

measurement errors, and 



Apj^^, Apz^ - 

filter outputs. 


The indicated along- 
components given in 

track is computed from the 
Equations (A. 52) and (A. 53) 

position error 
as 


ATD = AXj^ cost|j + AXg sinij; 

CA. 57) 

where 

is the nominal aircraft heading. 



For completeness it is desirable to represent the north and 
east displacements of Equations (A. 46) and (A. 47) in terms of 
various latitudes and longitudes. Using the primary facility 
as the origin of the coordinate system the various displace- 
ments are represented by 
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^N 


(A. 58) 



(A. 59) 

^2N 

= RpCtf)^ -tf>- ) , and 
e 52 s 

(A. 60) 

^2E 

= R^ cos(f)^ (X^ -X^ ) 
e ^S, St 

(A. 61) 


where 

((f) , A. } - aircraft latitude and longitude, 

3 - 3 

(cf) ,X ) - primary facility latitude and longitude, 

Si Sj 

and 

((f) ,X ) - secondary facility latitude and longitude. 
A. 3 MLS SIMULATION 

Figure A. 4 depicts the MLS simulation. For the current 
study it is assumed that on entering the MLS region that 
the navigation mode changes from RNAV to MLS without any 
blending. This implies that at the MLS boundary the navig- 
ation mode is switched directly and that the non-navigation 
elements of the covariance matrix are retained and the 
navigation elements are reinitialized as appropriate. The 
non-navigation states are those relating to the autopilots 
and aircraft dynamics. Hence, this subsection will not deal 
with the autopilot or aircraft dynamics as these have al- 
ready been discussed in the previous section. 

The primary measurements for the MLS consist of range, 
bearing and elevation angle. These measurements are refer- 
enced with respect to a coordinate system fixed at the run- 
way. Figure A. 5 indicates a schematic of the MLS measure- 
ments. Note that the bearing is with respect to the runway 
centerline rather than with respect to magnetic north as 
is the case for the VOR system. 

As with the RNAV system the MLS measurements are 
assumed to consist of a random noise and bias component. 
Again, for purposes of this study no correlated noise 
components are induced but can be easily added for subseq- 
uent studies. The measurement equations are given by 
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A, 4 MLS S iiTiul at. ion Block Diagram 
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Figure A. 5 MLS Antenna Coordinate System 



Pa 

* ^ 


(A. 62) 

^ = 

0 

a 

* ^ 

+ rig , and 

(A. 63) 

e = 
m 

£ 

a 

+ b 

e 

+ n , 
e ’ 

(A. 64) 


where (p ,0 ,e ) are the actual range, azimuth and elevation 

SL 3. 

angles, (b ,b , b ) are the bias components with their dy- 
namic behavior governed by 

^p = ^0 = ° 

and (rip > rig, n^) are the noise components. Linearization 
about the nominal yields 


Pm 

= Pn 

m 


(A. 66) 

0 

m 


+ A0 , and 
m’ 


(A. 67) 

e 

m 

= "n 

+ A £ , 
m’ 


(A. 68) 

Pa 

== Pn 

+ Ap , 


(A. 69) 

0 

a 


+ A0 , and 


(A. 70) 

e 

a 


+ Ae 


(A. 71) 

Substitution 

of these six 

equations into Equations 

(A. 62) 

through 

(A. 64) yields 



Ap 

m 

= Ap 

+ b + n , 
P P 


(A. 72) 

m 

= A0 

* '=8 * "8 ’ 

and 

(A. 73) 

Ae 

m 

= Ae 

* \ * ”e’ 

• 

(A. 74) 


The variation in the actual range, bearing and elevation 
with respect to the nominal, Ap , A0 , and Ac 
respectively are based on the actual aircraft cross-track 
and vertical deviations, CTD and Az. As with the RNAV 
simulation, the actual aircraft along-track deviation is 
not defined as 4D is considered beyond the scope of this 
study. The appropriate transformation from the actual 
aircraft deviations to the MLS measurement coordinate sys- 
tems results in the following relationships. 
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X y 

Ap = CTD ( ^ sin(0^ - 4 , ) - ^ cos(0p -li^) | + i Az (A. 75) 


A 0 = 


CTD 


/x 2 


cos ( 0 j^ - ip -e^) , and 


(A. 76) 


Ac = 


CTD 


Uxj;jXgl)^-XN^ + z2] Kx-Xg^)^+y^ 


(XR-Xgl)sinC9^-1,) * 


+ cosC0j^-^)j + 


COSSj^ 

■ — Az * 

n 72 2 2 


(A. 77) 


Here, x„ , y„ are the nominal aircraft lateral coordinates 

relative to the runway based MLS, z is the nominal air- 
craft altitude relative to the runway, 9j^ is the orient- 
ation of the runway centerline with respect to magnetic 
north, \jj is the nominal aircraft heading, Pj^, 0 j^ 

and e., are the nominal MLS measurements and x , is 
N el 

the displacement of the elevation antenna from the DME 
and azimuth antennas which are assumed to be located at 
the origin of the coordinate system. Note that the air- 
craft altitude is derived from the MLS measurements thereby 
eliminating the need for an altimeter model. 

To utilize the MLS measurements in the second order 
complementary filter requires a coordinate transformation 
from the polar p, 9 and e coordinate system to a 
rectangular runway referenced coordinate system. The trans- 
formation equations are given by 
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Ax 


R 


A/r = 


2 i/b -AC 
tan© Ax 




"R 


AB + 




AC 


2 )F 


■AC 


'R 


R 




A© 


AC, (A. 78) 


(A. 79) 


Az 


CxR-Xei)tanej^ 


>^R 


R 


fi 


Ax^ + 






R 


Ay 


R 


^R-^el^ "yR 




cos 


727 T 

^ ^yR 


Ae 

> 


where 


A = 1 + tan^Qj^-t- tan^Cj^-*- tan^^ tan^Sj^ , 

B - tan^£^ 


~ Xgi tan^Ej^-p^ , 


2 tane^ 

AA = 2 Ae + 

cos 0J^ COS 


2 tane 


■N 


cos 0J^ COS Ej^ 


AB = 


2 =^el 


2 

cos £ 


Ae, and 


N 


AC 


2 tane^ 


2 

cos E 


N 


Ae - 2pj^Ap 


(A. 80) 


(A. 81) 
(A. 82) 
CA.83) 

(A. 84) 


(A. 85) 


(A. 86) 
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For the MLS the second order complementary filter re- 
quires additional acceleration inputs. The accelerometers 
are assumed to be body-mounted hence a suitable attitude 
reference system is required for the appropriate coordinate 
transformation for suitable filter input. The linearized 
accelerometer error quantities are given by 


z = 

K = 


a + a il/ 

X y ^z 


a 

z '^y 


where 


a + a - a ilj 

y z X ^z 

a + a - a lb 

z X ^y y ^x 

• • • • • 

= ]h = a = a = a 

X y z 


= 0 




accelerometer errors 
misalignment errors, and 
accelerometer outputs. 


CA. 87) 
(A. 88) 
(A. 89) 
CA.90) 


The flight times of concern for this particular study are 
short compared to accelerometer and alignment drifts. Hence, 
these error quantities are assumed to be biases governed by 
the differential equations given by Equation. (A. 90). 

The MLS second order complementary filter operates in 
the coordinate system referenced at the runway. Hence, the 
accelerometer outputs must be rotated from the east/north 
system to the runway fixed system. The runway system is 
fixed hence the accelerations can be written as 

Xj^ — Xj^ cosOj^ Xg sin6|^, C-^* ^^) 

yj^ = Xj^ sinBj^ - x^ cos9j^, and (A. 92) 



where 6^^ is the runway orientation relative to magnetic 
north. 

Having resolved the accelerations into the runway co- 
ordinate system, the MLS second order complementary filter 
equations can be written as 
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(A. 94) 


^2 = Xj + Kj Cx^ - x^) H. Xp , (A. 95) 

^3 “ *'5 (='m ■ ='!>■ ■ CA.96) 

^4 ■= Xj + K3 (y^ - X4), (A. 97) 

^5 = Xj * K2 Cy^ - X4) * y^, (A. 98) 

^6 ' "5 (>'in ■ *4> . (A-S9) 

X? = Xj * Kj (Zjj - x^) , (A. 100) 

^8 ~ ^9 ^ ^2 ^ ^R> (A. 101) 

=^9 ' ''5 ■ ==7) > 

where 


/V /\ /s. /s /\ 

X = x^, y = X 4 , z = Xy, X = X 2 » y = Xg , 

K 2 , Kj, K^- filter constants, 

X , y , z - MLS measurements, and 
m m m 

Xr, yR, Zr- given by Equations (A. 91) through (A. 93). 

Several coordinate transformations are required to 
resolve the runway referenced filter outputs into values 
appropriate for input into the steering command comput- 
ations. The first coordinate transformation resolves the 
estimated velocity into a heading error estimate given by 




sinif; ^ cosifj " 
Vn "n "e 


CA.103) 


where Vj^ is the nominal aircraft velocity and the north 

and east velocity error estimates are given in terms of the 
runway reference velocity error estimates as 
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Xn = Xj^ COSQ^ + Yj^ Sin0j^ , 


and 


E 


7r cos9j^ . 


CA. 104) 
(A. 105) 


where 9^ is the orientation o£ the runway relative to magnetic 
north. The indicated cross-track deviation CTD is given by 


CTD = - Xj^ sinij; + x^ cosi|) 


CA. 106) 


where tp is the nominal aircraft heading and the north and 
east position error estimates are given in terms of the 
filter output as 


^ ^R " T"r sinOj^, and 


Xg = Xj^ sinOj^ 


CA. 107) 
CA. 108) 


7r cosOj^ . 

The indicated altitude error is simply given by 

CA.109) 

is the filter estimate of the altitude error. 


>'c ' 


where 


'R 


A. 4 SIMULATION STATES 


The states of the VOR/DME RNAV system simulation are 
^RNAV ■tCTDjijj, z, h^, (p , (p^, h^, h^, h^^ , Xj^, 

^2N’ ^E’ ^2E’ ^^E’ *^^E ’ ’^x’ '^y’ 

i|^z> £y, £^, a^, Uy, b^} 

CA. 110) 

With the differential equations expressed in the form 

^RNAV " ^RNAV^ ^ CA.lll) 

it can be observed that the F matrix elements are made up 
of the equations presented in Section A. 2. For the DME/DME 
RNAV simulation the states are identical as the VOR/DME with 
the exception that bg is replaced by b 
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For this study it is assumed that the state transition 
matrix can be expressed by the following approximation 

<)> = I + ^RNAv'^^ (A. 112) 

where I is the identity matrix, F is the matrix des- 
cribing the dynamics of the system and At is the integration 
step size.' 

The formulation for the MLS simulation is the same with 
the exception that the MLS states are given by 


ms 


= {CTD, 


z, 

>'1 

> <i> J 



^2 ’ 

^1’ ^2’ ^3’ 


^5’ 

^6’ 

X? 

X 

00 

^9 > 

’^x’ 

ij; , 

, '^x’ “y’ 









(A. 113) 

matrix 

associated 

with 

the 

MLS 

is 

given by the 


The F..y o 
MLS 

equations presented in Section A. 3, 


The measurement noise matrix for the VOR/DME RNAV 
system is 


Q 


1 



0 


0 


0 



(A. 114) 


and for the DME/DME RNAV system as 


o 


2 

P 


Qz 


0 


0 

2 

a 

P 


0 0 



(A. 115) 
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where 


Oq - VOR bearing error variance. 

Op - DME range error variance, and 

a - altimeter error variance, 
n 

Similarly, the MLS measurement matrix is 



where 

o ^ - MLS range error variance, 

P 

- MLS azimuth error variance, and 
u 

o - MLS elevation angle error variance, 
e 


(A. 116) 
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APPENDIX B 

RNAV/MLS SIMULATION GEOMETRIES 


The purpose of this appendi:x is to. present the details 
of the geometries for the simulation experiments. The geo- 
metries are shown in Figures B.l, B.2 and B.3 for the DME/ 
DME analysis and Figures B.4, B.5 and B.6 for the VOR/DME 
analysis. Tables B.l, B.2, and B.3 contain the waypoint 
and station data for DME/DME and Tables B.4, B.5 and B.6 
contain the waypoint and station data for VOR/DME. The 
figures are duplicates of those presented in Section III 
of the report and are merely included here for reference. 

IN all cases the runway length is assumed to be 3050 meters. 

Table B.l 


Straight- In Approach for DME/DME Navigation Mode 



LAT 

LON 

ALTITUDE 

Waypoint WP.| 

39°00'00" 

73°27'11.3" 

2073 meters 

Waypoint WP 2 

39°00'00" 

73°30'24.3" 

2073 meters 

Waypoint WP^ 

39099 '00" 

73®57'57.4" 

0 meters 

Facility a^ 

38°56'10.4" 

73°22'23.3" 

0 meters 

Facility a 2 

38°56'10.4" 

73°46'09.2" 

0 meters 

Facility b^ 

39°07'04.3" 

73°43'09.2" 

0 meters 

Facility 62 

38°52'55.7" 

73°43'22.7" 

0 meters 
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Figure B.l Straight- In Approach for DME/DME 
Navigation Mode 
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Figure B.3 180° Turn to Final Approach for 

DME/DME Navigation Mode 
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Straight 
and Level 



3° Descent 


Q 

VOR/d:i: Facility 


Figure B.4 Straight-In Approach for VOR/DME Navigation 
Mode 




Figure B.5 90® Turn to Final VOR/DME Navigation Mode 




Figure B. 6 180° Turn to Final for VOR/DME Navigation 

Mode 
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Table B . 2 


'90° Turn- to- Final Approach for DME/DME 
Navigation Mode 



LAT 

LON 

ALTITUDE 

Waypoint WP^ 

38°49'24.4" 

73°52'11.1" 

1165 meters 

Waypoint WP 2 

38°51 '54.4" 

73°52'11 .1" 

1165 meters 

Waypoint WP^ 

38°58'46.0" 

73°52'11.1" 

499 meters 

Waypoint WP^ 

39°00'00" 

73°53'46.3" 

311 meters 

Waypoint WPg 

39°00'00" 

73°57'53.4" 

0 meters 

Facility a^ 

39°04'08.7" 

73°57'07.5" 

0 meters 

Facility a^ 

38°45'40.1" 

73°57'07.5" 

0 meters 

Facility 

38°58'44.0" 

74°04'04.2" 

0 meters 

Facility b 2 

38°58'^4.0" 

73°40'18.0" 

0 meters 


Table B.3 

180° Turn-to-Final Approach for DME/DME 
Navigation Mode 



LAT 

LON 


Waypoint WP^ 

38° 56 '42.6" 

74°00'50.3" 

853 meters 

Waypoint WP 2 

38°56'42.6" 

73°57'37.4" 

853 meters 

Waypoint WP^ 

38°56'42.6" 

73°53'45.9" 

526 meters 

Waypoint WP^ 

39°00'00" 

73°53'45.9" 

311 meters 

Waypoint WPg 

39°00'00" 

73°57'53.4" 

0 meters 

Facility a^ 

38°52'53.0" 

74°05'38.0" 

0 meters 

Facility a 2 

38°05'56.9" 

73°58'41.8" 

0 meters 

Facility 

39°05'56.9" 

73°58'41.8" 

0 meters 

Facility b 2 

38°47'28.3" 

73°58'41.8" 

0 meters 
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Table B.4 

Straight- In Approach for VOR/DME Navigation Mode 



LAT 

LON 

ALTITUDE 

Waypoints WP^ , 
Facility A 
Facility B 

WP 2 . WP 3 
38°40'00" 
39°00'00" 

Same as T 
73°34'15.9" 
73°08'31 .8" 

Jble B.l 

0 meters 
0 meters 


Table B.5 

90° Turn- to- Final Approach for VOR/DME 
Navigation Mode 



LAT 

LON 

ALTITUDE 

Waypoints WP-j , 
WP4, 

Facility A 
Facility B 

WP2. WP 3 , 

WP5 

38° 54' 54.9" 
38°34'54.9" 

Same as T; 

74°41'21.0" 

73°52'11.1" 

Jble B.2 

0 meters 
0 meters 


Table B.6 

180° Turn- to- Final Approach for VOR/DME 
Navigation Mode 



LAT 

LON 

ALTITUDE 

1 Waypoints WP^ , WP 2 , WP 3 , 

1 

Same as Table B.3 

WP4. 1 

■">5 



Facility A 

39°33’49.6" 

73°53’45.9" 

0 meters 

Facility B 

38°56'42.6" 

74°19.28.9" 

0 meters 
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